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Filling  N-sided  Holes  Using 
Combined  Subdivision  Schemes 


Adi  Levin 


Abstract.  A new  method  is  presented  for  calculating  TV-sided  surface 
patches.  The  algorithm  generates  a subdivision  surface  which  satisfies 
arbitrary  C 1 boundary  conditions.  The  proposed  subdivision  scheme  is 
based  on  a Catmull-Clark  type  subdivision  scheme  that  operates  in  the 
surface  interior.  Near  the  boundary  we  introduce  new  subdivision  rules 
that  involve  the  given  boundary  conditions.  The  generated  subdivision 
surface  is  C2-continuous  except  at  one  extraordinary  point.  In  the  neigh- 
borhood of  this  point  the  surface  curvature  is  bounded. 


§1.  Background 

The  problem  of  constructing  TV- sided  surface  patches  occurs  frequently  in 
computer-aided  geometric  design.  The  TV-sided  patch  is  required  to  connect 
smoothly  to  given  surfaces  surrounding  a polygonal  hole,  as  shown  in  Fig.  1. 

Referring  to  [10,25,26],  TV-sided  patches  can  be  generated  basically  in 
two  ways.  Either  the  polygonal  domain,  which  is  to  be  mapped  into  3D,  is 
subdivided  in  the  parametric  plane,  or  one  uniform  equation  is  used  to  repre- 
sent the  entire  patch.  In  the  former  case,  triangular  or  rectangular  elements 
are  put  together  [2,6,12,20,23]  or  recursive  subdivision  methods  are  applied 
[5,8,24].  In  the  latter  case,  either  the  known  control-point  based  methods  are 
generalized,  or  a weighted  sum  of  3D  interpolants  gives  the  surface  equation 
[1,3,4,22]. 

The  method  presented  in  this  paper  is  a recursive  subdivision  scheme 
specially  designed  to  consider  arbitrary  boundary  conditions.  Subdivision 
schemes  provide  efficient  algorithms  for  the  design,  representation  and  pro- 
cessing of  smooth  surfaces  of  arbitrary  topological  type.  Their  simplicity  and 
their  multiresolution  structure  make  them  attractive  for  applications  in  3D 
surface  modeling,  and  in  computer  graphics  [7,9,11,13,19,27,28]. 

The  subdivision  scheme  presented  in  this  paper  falls  into  the  category 
of  combined  subdivision  schemes  [14,15,17,18],  where  the  underlying  surface  is 
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Fig.  1.  A 5-sided  and  a 3-sided  surface  patch. 

represented  not  only  by  a control  net,  but  also  by  the  given  boundary  con- 
ditions. The  scheme  repeatedly  applies  a subdivision  operator  to  the  control 
net,  which  becomes  more  and  more  dense.  In  the  limit,  the  vertices  of  the 
control  net  converge  to  a smooth  surface.  Samples  of  the  boundary  conditions 
participate  in  every  iteration  of  the  subdivision,  and  as  a result  the  limit  sur- 
face satisfies  the  given  conditions,  regardless  of  their  representation.  Thus, 
our  scheme  performs  so-called  transfinite  interpolation. 

The  motivation  behind  the  specific  subdivision  rules,  and  the  smoothness 
analysis  of  the  scheme  are  presented  in  [16].  In  the  following  sections,  we 
describe  Catmull-Clark’s  scheme,  and  we  present  the  details  of  our  scheme. 

§2.  Catmull-Clark  Subdivision 

A net  E = (V,E)  consists  of  a set  of  vertices  V and  the  topological  information 
of  the  net  E,  in  terms  of  edges  and  faces.  A net  is  closed  when  each  edge  is 
shared  by  exactly  two  faces. 

Camull-Clark’s  subdivision  scheme  is  defined  over  closed  nets  of  arbitrary 
topology,  as  an  extension  of  the  tensor  product  bi-cubic  B-spline  subdivision 
scheme  [5,8].  Variants  of  the  original  scheme  were  analyzed  by  Ball  and  Storry 
[24].  Our  algorithm  employs  a variant  of  Catmull-Clark’s  scheme  due  to  Sabin 
[21],  which  generates  limit  surfaces  that  are  C2-continuous  everywhere  except 
at  a finite  number  of  irregular  points.  In  the  neighborhood  of  those  points 
the  surface  curvature  is  bounded.  The  irregular  points  come  from  vertices  of 
the  original  control  net  that  have  valency  other  than  4,  and  from  faces  of  the 
original  control  net  that  are  not  quadrilateral. 

Given  a net  E,  the  vertices  V'  of  the  new  net  E'  = (V',  E')  are  calculated 
by  applying  the  following  rules  on  E (see  Fig.  2): 

1)  For  each  old  face  /,  make  a new  face-vertex  v(f)  as  the  weighted 
average  of  the  old  vertices  of  /,  with  weights  Wm  that  depend  on  the  valency 
m of  each  vertex. 

2)  For  each  old  edge  e,  make  a new  edge-vertex  v(e)  as  the  weighted 
average  of  the  old  vertices  of  e and  the  new  face  vertices  associated  with  the 
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Fig.  2.  Catmull-Clark’s  scheme. 


two  faces  originally  sharing  e.  The  weights  Wm  (which  are  the  same  as  the 
weights  used  in  rule  1)  depend  on  the  valency  m of  each  vertex. 

3)  For  each  old  vertex  v,  make  a new  vertex-vertex  v(v)  at  the  point  given 
by  the  following  linear  combination,  whose  coefficients  atm,pm,'ym  depend  on 
the  valency  mofu: 

am ■ (the  centroid  of  the  new  edge  vertices  of  the  edges  meeting  at  v)  + 
pm-  (the  centroid  of  the  new  face  vertices  of  the  faces  sharing  those  edges)  + 
7 m-v. 

The  topology  E'  of  the  new  net  is  calculated  by  the  following  rule:  for 
each  old  face  / and  for  each  vertex  v of  /,  make  a new  quadrilateral  face 
whose  edges  join  v(f)  and  v(v)  to  the  edge  vertices  of  the  edges  of  / sharing 
v (see  Fig.  2). 

We  present  the  procedure  for  calculating  the  weights  mentioned  above, 
as  formulated  by  Sabin  in  [21]:  Let  m > 2 denote  a vertex  valency.  Let 
k :=  cos(7T /m).  Let  x be  the  unique  real  root  of 

Z3  + (4 k2  - 3)x  -2k  = 0, 

satisfying  x > 1.  Then 

Wm  = x2  + 2 kx  - 3,  am  = 1, 

kx  2k 2 — 1 _ 

7m  ‘ o / 7 j Pm  7 m- 

xz(kx  + 1) 

The  weights  Wm  and  7m  for  m = 3, . . . , 7 are  given  by 


W3  = 1.23606797749979 . . . 
W4  = 1.00000000000000000 
W5  = 0.71850240323974... 
W6  = 0.52233339335931 . . . 
W7  = 0.39184256502794 . . . 


73  = 0.06524758424985 . . . 

74  = 0.25000000000000000 

75  = 0.40198344690335 . . . 

76  = 0.52342327689253 . . . 

77  = 0.61703187134796 . . . 


Remark:  The  original  paper  by  Sabin  [21]  contains  a mistake:  the  formulas 
for  the  parameters  a,p  and  7 that  appear  in  §4  there  are  p :=  1,  7 :=  —a. 
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Fig.  3.  The  input  data  (left)  and  the  initial  control  net  (right). 

§3.  The  Boundary  Conditions 

The  input  to  our  scheme  consists  of  N smooth  curves  given  in  a parametric 
representation  Cj  : [0,2]  — » IR3  over  the  parameter  interval  [0,2],  and  corre- 
sponding cross-boundary  derivative  functions  dj  : [0,2]  — » 1R3  (see  Fig.  3). 
We  say  that  the  boundary  conditions  are  C°-compatible  at  the  j-th  corner  if 


cj(2)  — cj+i(0). 


We  say  that  the  boundary  conditions  are  C^-compatible  if 

dj(  0)  = -<$_,(  2), 

djW  = c'j+1(0). 

We  say  that  the  boundary  conditions  are  C2-compatible  if  the  curves  cj  have 
Holder  continuous  second  derivatives,  the  functions  dj  have  Holder  continuous 
derivatives,  and  the  following  twist  compatibility  condition  is  satisfied: 

d'j(2)  = -d’j+1(0).  (1) 

The  requirement  of  Holder  continuity  is  used  in  [16]  for  the  proof  of  C 2- 
continuity  in  case  the  boundary  conditions  are  C2-compatible. 

§4.  The  Algorithm 

In  this  section  we  describe  our  algorithm  for  the  design  of  an  IV-sided  patch 
satisfying  the  boundary  conditions  described  in  §3.  The  key  ingredients  of  the 
algorithm  are  two  formulas  for  calculating  the  boundary  vertices  of  the  net. 
These  formulas  are  given  in  §4.3  and  §4.4. 
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Fig.  4.  The  stencils  for  the  smooth  boundary  rule  (left)  and  the  corner  rule 
(right). 


4.1.  Constructing  an  initial  control  net 

The  algorithm  starts  by  constructing  an  initial  control  net  whose  faces  are  all 
quadrilateral  with  2 N boundary  vertices  and  one  middle  vertex,  as  shown  in 
Fig.  3.  The  boundary  vertices  are  placed  at  the  parameter  values  0, 1, 2 on 
the  given  curves.  The  middle  vertex  can  be  arbitrarily  chosen  by  the  designer, 
and  controls  the  shape  of  the  resulting  surface. 

4.2.  A single  iteration  of  subdivision 

We  denote  by  n the  iteration  number,  where  n = 0 corresponds  to  the  first 
iteration.  In  the  n-th  iteration  we  perform  three  steps:  First,  we  relocate  the 
boundary  vertices  according  to  the  rules  given  below  in  §4.3  - §4.4.  Then, 
we  apply  Sabin’s  variant  of  Catmull-Clark’s  scheme  to  calculate  the  new  net 
topology  and  the  position  of  the  new  internal  vertices.  For  the  purpose  of 
choosing  appropriate  weights  in  the  averaging  process,  we  consider  the  bound- 
ary vertices  as  if  they  all  have  valency  4.  This  makes  up  for  the  fact  that  the 
net  is  not  closed.  In  the  third  and  final  step,  we  sample  the  boundary  ver- 
tices from  the  given  curves  at  uniformly  spaced  parameter  values  with  interval 
length  2-(n+1). 


4.3.  A smooth  boundary  rule 

Let  v denote  a boundary  vertex  corresponding  to  the  parameter  0 < u < 2 on 
the  curve  cj.  Let  w denote  the  unique  internal  vertex  which  shares  an  edge 
with  v (see  Fig.  4 (left)).  In  the  first  step  of  the  n-th  iteration  we  calculate 
the  position  of  the  v by  the  formula 

v = 2c, (it)  - ^ {cj  {u  + 2"")  + Cj  ( u - 2~n)) 

- 2-"i  (dj  (u  + 2~n)  + dj  (it  - 2-"))  + 2-"^(it). 


226 


Adi  Levin 


4.4.  A corner  rule 


Let  v denote  a boundary  vertex  corresponding  to  the  point  c_,_i(2)  = Cj(0). 
Let  w be  the  unique  internal  vertex  sharing  a face  with  v (see  Fig.  4 (right)). 
In  the  first  step  of  the  n-th  iteration  we  calculate  the  position  of  v by  the 
formula 


5c.(o)  + lw  - (cj(2~n)  + c_,_i(2  - 2-))  + \cj{2l~n) 


,29 


1 


+ -Cj-1(2  - 21_n)  + 2-"-  (dj (0)  + dj_i(2))  - -d,{ 2~B) 


48 


12 


• 2-»idj_1(2  - 2-)  - 2-»l  (dj(21-n)  + dj -i(2  - 21""))  . 


§5.  Properties  of  our  Scheme 

In  [16]  we  prove  that  the  vertices  generated  by  the  above  procedure  converge  to 
a surface  which  is  C2-continuous  almost  everywhere,  provided  that  the  bound- 
ary conditions  are  G2-compatible  (as  defined  in  §3).  The  only  point  where  the 
surface  is  not  C2-continuous  is  a middle-point  (corresponding  to  the  middle 
vertex,  which  has  valency  N),  where  the  surface  is  only  G^continuous.  In  the 
neighborhood  of  this  extraordinary  point,  the  surface  curvature  is  bounded. 

The  limit  surface  interpolates  the  given  curves,  for  C°-compatible  bound- 
ary conditions.  For  G^cornpatible  boundary  conditions,  the  tangent  plane  of 
the  limit  surface  at  the  point  Cj(u)  is  spanned  by  the  vectors  c'-(u)  and  dj(u), 
thus  the  surface  satisfies  G^boundary  conditions.  Furthermore,  due  to  the 
locality  of  this  scheme,  the  limit  surface  is  G2  near  the  boundaries  except  at 
points  where  the  G2-compatibility  condition  is  not  satisfied. 

The  surfaces  in  Fig.  5 demonstrate  that  the  limit  surface  behaves  mod- 
erately even  in  the  presence  of  wavy  boundary  conditions.  The  limit  surfaces 
are  G2-continuous  near  the  boundary  except  at  corners  where  the  twist  com- 
patibility condition  (1)  is  not  satisfied. 


Fig.  5.  A 5 -sided  and  a 3-sided  surface  patch  with  wavy  boundary  curves. 
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